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In our work we have reconsidered the old problem of diffusion at the boundary of 
ultrametric tree from a "number theoretic" point of view. Namely, we use the modu- 
lar functions (in particular, the Dedekind 77-function) to construct the "continuous" 
analog of the Cayley tree isometrically embedded in the Poincare upper half-plane. 
Later we work with this continuous Cayley tree as with a standard function of a com- 
plex variable. In the frameworks of our approach the results of Ogielsky and Stein on 
dynamics on ultrametric spaces are reproduced semi-analytically/semi- numerically. 
The speculation on the new "geometrical" interpretation of replica n — > limit is 
proposed. 
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I. INTRODUCTION 



Let us begin with an obvious and almost tautological statement: any regular tree has 
an ultrametric structure. Recall that ultrametricity of the space H, implies the "strong 
triangle inequality" meaning that the distances tab, Tbc and vac between any three points 
A, B,C in the space Ti satisfy the condition tag < max \jab-, r Bc} ( see , f° r example, jl|). 
Appearing in mathematical literature in connection with p-adic analysis (see, for review 



Appef 

HB) 



ultrametric spaces became very popular in physical community because of their 
importance for spin glasses — see, for review f4|, and references therein. 

The famous replica symmetry breaking (RSB) scheme is ultimately connected to the 
ultrametric structure of the phase space of many disordered systems possessing the spin- 
glass behavior. Since invention of RSB, many authors, both physicists and mathematicians, 
attempted to adapt the p-adic analysis for physical needs mainly trying to elucidate and 
justify the replica n — > limit. However to our point of view the interpretation of spin- 
glass problems in terms of p-adic language does not converge well. As an exception one 
has to mention recent interesting contributions to that subject 0, 0, H| , where the appli- 
cation of p-adic Fourier transform allows significantly simplify the solution to the problem 
of diffusion on ultrametric tree. One can hope that the continuation and generalization of 
these works would release deeper penetration of p-adic analysis into physics of disordered 
systems with ultrametric phase spaces, such as spin glasses, neural networks and disordered 
heteropolymers (see |H E| for the detailed list of corresponding references) . 

In what follows we shall always keep in mind the p + 1-branching Cayley tree as an 
example of an ultrametric space. Despite extremely simple topological structure of a tree, 
one cannot operate in this space as in usual space with an Euclidean metric because the 
number of degrees of freedom for the p + 1-branching Cayley tree grows exponentially with 
the size of a tree. For some problems, such as, for example, the branching process on Cayley 
trees and tree-like graphs it is not sufficient to deal only with the "distance" measured in 
number of generations between two points on the graph, but it is ultimately necessary to 
know the absolute values of coordinates of points on the Cayley tree. The main difficulty 
concerns the encoding of the Cayley tree vertices. This problem becomes very cumbersome 
because for the tree we do not have any transparent and convenient "coordinate system", 
like a D-dimensional grid in a D-dimensional Euclidean space. One of possible ways to 
resolve the addressed problem consists in using the p-adic analysis 11 1 . The vertices of 



the graph C, i.e. of the p + 1-branching Cayley tree, admit the natural parameterizations 
by the p-adic numbers what enables to develop the whole machinery like p-adic Fourier 
transform etc. This way has been exploited in papers 0,0, H[. 

Another possibility, described in the present paper consists in the following. Instead of 
working with the ultrametric discrete graph C, one can embed this graph in the metric 
space 7i preserving all ultrametric properties of C. Let us mention that any image of a 
regular Cayley tree refers indirectly to the isometric embedding of a tree into the complex 
plane. Indeed, the Cayley tree is usually drawn in the picture in such a way that any new 
generation of vertices (counted from the root point) is smaller than the previous generation 
in geometric progression. Taking advantage of the isometric embedding, one can: 



- naturally parameterize the vertices of the Cayley graph C by ordinary complex variable 
z = x + iy in the complex plane z without using any ingredients of a p-adic analysis; 
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- construct a continuous analog of a Cayley tree, i.e. a continuous space TL with ultra- 
metric properties borrowed from the initial Cayley tree. 

The paper is organized as follows. In the section |H] we describe the way of isometric 
embedding of the regular 3-branching Cayley tree C into the complex upper half-plane 
z = x + iy as well as we discuss the ultrametric properties of the resulting space TL. The 
problem of diffusion on the boundary of ultrametric tree C embedded into the space TL is 
considered in the section UTT1 The results are summarized in the conclusion, where we express 
also some conjectures concerning the possible geometrical interpretation of replica n — > 
limit. 



II. ULTRAMETRIC STRUCTURE OF ISOMETRIC CAYLEY TREES 

It is well known that any regular Cayley tree, as an exponentially growing structure, 
cannot be isometrically embedded in an Euclidean plane. Recall that the embedding of a 
Cayley tree C into the metric space is called "isometric" if C covers that space, preserving all 
angles and distances. For example, the rectangular lattice isometrically covers the Euclidean 
plane S = {x, y} with the flat metric ds 2 = dx 2 + dy 2 . In the same way the Cayley tree C 
isometrically covers the surface of constant negative curvature (the Lobachevsky plane) TL. 
One of possible representations of TL, known as a Poincare model, is the upper half-plane 
Im z > of the complex plane z = x + iy endowed with the metric ds 2 = dx ^ 2 dy of constant 
negative curvature. 

In this Section we are aimed to construct the "continuous" analog of the standard 3- 
branching Cayley tree by means of modular functions and analyse the structure of the 
barriers separating the neighboring valleys. Due to specific number-theoretic properties of 
modular functions these barriers are ultrametrically organized. 

A. Continuous analog of isometric Cayley tree 

To be precise, let us begin with the explicit description of the standard recursive con- 
struction which allows encoding of all vertices of the 3-branching Cayley tree C isometrically 
covering the surface of the constant negative curvature TL = {z\ Imz > 0}. Recall that the 
3-branching Cayley tree is the Cayley graph of the group A which has the free product 
structure: A ~ Z2 <S> ^2 ® ^2 (where Z2 is the cyclic group of 2nd order). The matrix rep- 
resentation of the generators hi,Ji2, of the group A is well known (see, for example jl2j): 

ft -(ri)' A *=(^)-''-Ut) « 

For our purposes it is convenient to take a framing consisting of the composition of the 
standard fractional-linear transform and the complex conjugacy. Namely, denoting by z the 
complex conjugate of z, we consider the following action in 7i: 

(a b \ az + b 

c dj Z ^ cz + d ^ ' 
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We raise the Cayley tree C in 7i as follows. Take the set of generators {hi,h2,hs} of the 
group A. Choose the point (x , iy ) = (0, i) as the root of the tree. 

Any vertex in the generation n from the root point of the tree C is associated with 

n 

an element M n = j^J h ak of the group A (where G {1,2,3} for any k). All vertices 

k=l 

are parameterized by the complex coordinates z n = M n ^(— l) n i^j in 7i. For example, the 

coordinates of the points Za and in figffl we can get the following way. Multiply the 
generators hi along the trajectory in the backway order, i.e. from the points za and zb to 
the root point i. Hence we arrive at 

M { 2 A) = I I = //, // 



Mi B) 





Using (J2J and the matrices and M^ B \ we get 

0x ^-73 2 

Z A = — 7= — = 7= + ■= 

(-1) x (-i) + 2^ % 
ZB ~ 2v/3x (-0 + 1 ~ ~^3" + 13 

To establish the connection with the forthcoming constructions, we shall make the simple 
linear transform oi z : z ^ + \. The coordinates of the points A and B are shown in 
FigQ after such transform. 

However such recursive construction has, to our point of view, the crucial lack: one cannot 
automatically encode all vertices of the tree introducing a sort of a " coordinate system" as, 
say, in the Euclidean plane. Actually, we know that any vertex of a rectangular grid in the 
Euclidean plane has the coordinates (am, bn), where {m,n} G Z and (a, b) are the length 
and the width of the elementary cell of the lattice. Our desire is to construct somewhat 
similar for homogeneous trees. To be more specific, we are aimed to find an analytic function 
f(z) defined in Im z > whose all zeros give all coordinates of the vertices of the 3-branching 
homogeneous Cayley tree C isometrically embedded into H.(z\lmz > 0). In the rest of this 
Section we describe the construction of the corresponding function f(z\lm z > 0) and discuss 
its properties in connection with the geometry of ultrametric spaces. The forthcoming 
construction is based on the properties of the Dedekind 77-function. Recall the standard 
definition of the rj(z) -function (see, for instance 

00 

v (z) = e wiz/12 Y[{1 - e 2wikz ); lmz>0 (3) 

fc=0 

It is well known 0] that the Dedekind ^-function is connected to the elliptic Jacobi $- 
functions by the following relation: 

$[(0,e mz ) = V 3 (z) (4) 



? 
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FIG. 1: Left: 3-branching Cayley tree isometrically embedded in Poincare hyperbolic upper half- 
plane H. Right: Density plot of the function f(z) (see the text) in the rectangle {0 < Rez < 
1, 0.01 < Imz < 1.4}. 



where 

(11 p niz \ JUL 
^(0, e niz ) = U J ' = 2e mz/4 ^(-l) n (2n + l) e ™( n+l > (5) 

u=0 n= o 

Now we are in position to formulate the central assertion of the paper. Define the function 
f(z) as follows: 

f(z) = C- 1 \r,(z)\Qmz) 1 '* 
where C is the normalization constant, 

1/4 



(6) 



C 



I = 0.77230184... 



(7) 



The normalization constant C is chosen to set the maximal value of the function f(z) to 1: 
< f(z) < 1 for any z in the upper half-plane Imz > 0. It can be proved that the function 
f(z) has the following remarkable properties: 

1. The function f(z) — 1 has zeros at the points z ccn (and only at these points) - the 
centers of zero-angled triangles tesselating the Poincare hyperbolic upper half-plane 
H{z\Imz > 0); 

2. The function f(z) has local maxima at all the points z cen and only at these points. 

All the solutions of the equation f(z) — 1 = define all the coordinates of the 3-branching 
Cayley tree isometrically embedded into the upper half-plane TC(z\lmz > 0). The corre- 
sponding density plot of the function f(z) in the region {0 < Re z < 1, 0.04 < Im z < 1.4} is 
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shown in figO It is noteworthy that the function Z(z) = [C f(z)}~ 2 is the partition function 
of a free bosonic field on a torus 1.14.] • The function Z(z) has been extensively studied in 
the conformal field theory. In particular it can be easily checked that the function Z(z) is 
invariant with respect to the action of the modular group PSL(2, Z), namely: 



Z(z) =Z(z + l] 

1\ (8) 

z 



Z{z) = Z 



Later on we shall address to these self-dual properties. The proof of our assertion © implies 
the check that the function f(z) is invariant with respect to the conformal transform z^(z) 
of the Poincare half-plane to itself 

2<1, W = tFt < 9 > 

ZZq — 1 

where z is the coordinate of any center of zero-angled triangle in the hyperbolic Poincare 
upper half-plane obtained by successive transformations from the initial one. Hence, it is 
necessary and sufficient to show that the values of the function f(z) at nearest centers of 
circular triangles are equal and reach the maximal value 1: 

Then taking z = { (j + ^"^r ) > ( — I + > (l + ^^ff) } an< ^ performing the conformal 
transform Q we move the centers of the first generation of zero-angled triangles to the new 
centers (second generation) located at z^ 1 '. Now we can repeat recursively the construction, 
i.e. find the new coordinates of the centers of the third generation of zero-angled triangles, 
and compute the function f(z) at these points, then we perform the conformal 

transform 

2 (3)( 2 (2)) 

and so on... Let us use now the following well-known properties of Jacobi 
elliptic ^-functions [13] 

pz + r 



7] ( = uy/rz + s 7](z) (11) 

\qz + s J 

where Im z > 0, {p, q, r, s} G Z; ps — qr = 1 and u) is some 24-power root of unity, which 
depends on the coefficients p, q, r, s but does not depend on z. Using (fTTj) we can write: 

r}(z + k) = r]{z)] fc = ±l,±2,±3,... 

v( 1 , + ^)^ 1/4 = vU + ^); A>o 



2^2^/ 1 \2 ^ 2 

Now we can easily check the relations (JTUJ). In turn, (jHJ) can be regarded as a particular 
case of eq.([ll|). The fact that the function f(z) has local maxima at the points z cen can be 
verified straightforwardly by differentiating f(z) with respect to z at the points z cen . The 
3D relief of the function f(z) is shown in figure El The tree-like structure of hills separated 
by the valleys is clearly seen. 
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FIG. 2: Relief of the function f(z) in the rectangle {0 < Rez < 1, 0.04 <Imz< 1.4}. 



Let us note the important duality relation for the Dedekind function which establishes 
the mapping y «-> ^ (where y = Im z; y > 0). For all {p, g, r, s} G 2" such that ps — qr = 1 
the following equation is satisfied: 



ly 



s 



q 2 y 



1 



(13) 



where by ~ > and 1 1 1 we denote the fractional parts of the corresponding quotients. The 

relation (fl~3*|) shall be used in the forthcoming sections for practical purposes of numerical 
computations. In the next chapter we discuss the ultrametric organization of corresponding 
barriers. We regard the 3D relief constructed on the basis of the function f(z) (see figEJ) as 
the continuous analog of the Cayley tree isometrically embedded in Ti,(z\lmz > 0). 



B. Ultrametric structure of barriers in the standard RSB scheme and in the 

tree like metric space 

The replica symmetric breaking (RSB) scheme has appeared in the spin-glass theory 
as a self-consistent approach free of many lacks of the replica symmetric consideration. Later 
on it has been realized that the RSB-structure naturally appears in the problem of diffusion 
on the boundary of the Cayley tree 0, 0], where the neighboring sites are separated by 
the barriers hierarchically organized according their ultrametric distances on the Cayley 
tree. For example the neighboring points A, B and B, C at the boundary of the 3-branching 
Cayley tree — see figO] — are separated by the barriers depending on the ultrametric distance 
between points A, B and B, C on the tree. 



8 



Recall that two points, say, B and C in figE] are separated by the potential barrier 
depending on the number of the Cayley tree generations from the points B and C up to 
their common "parent branch". Let us label the states at the L's generation of the 3- 
branching Cayley tree by the integer number k (1 < k < 2 L ). If we represent k by the 
binary sequence, then the ultrametric distance between points B and C shall coincide with 
the highest distinct rank in the binary writing of ks = 4 and kc = 5. 




FIG. 3: The 3-branching Cayley tree. The ultrametric distance between points A, B and B, C are 
shown by bold lines; n is the number of states in each generation of the tree, the maximal number 
of generations L = 4. 

Denote by q m the Boltzmann weights at the boundary of 3-branching Cayley associated 
to the barriers of ultrametric height m. All the barriers can be encoded in the 2 L x 2 L RSB 
matrix Q(9o,9i, ■■■,qi J ) having the structure shown in eq- lfPIj) for L = 3. Keeping in mind 
the connection to the spin glasses, the total number of states 2 L is associated to the total 
number of replicas, n, in the standard RSB scheme. 



Q 
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92 92 
92 92 


90 9i 

91 9o 



(14) 



The values p m = \q m ( m — 1 ; •••) L) define the jumping probabilities from some vertex ki 
to another vertex &2 at the boundary of 3-branching Cayley tree separated by the barrier 
of ultrametric height m, while the value po is assigned for the probabilities to stay in a 
given vertex. In each line of the matrix Q the sum of probabilities is equal to 1. In most 
physically important situations the Boltzmann weights q m depend on m either exponentially 
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or polynomially: 



Qn 



const m a 
const e~^ m 



(15) 



Later on we pay attention only to the exponential dependence of q m on m, i.e. to the 
case when q m = const e _/3m (where (3 stands for the inverse temperature). The probability 
Po to stay in a given vertex can be explicitly computed from the conservation condition 
Ylm=o P™, = 1 and for q m = const e~ 13171 it reads: 



Pa 



n ^ 



-Pi 



2(1 -2 L e~P L ) 
n(2 - el 3 ) 



(16) 



In order to understand the connection between the ultrametric structure of the barriers 
in case of 3-branching Cayley tree and its continuous analog isometrically embedded in 
TC(z\lmz > 0), let us begin with the following observation. We can consider the function 
v(x\y) = — ha.f(x\y) (where x = Rez) as the potential relief at the boundary of the Cayley 
tree cut at the distance y = Imz from the real axis. The typical shapes of the function 
v(x\y) for few fixed values y = 0.04, 0.008, 0.0016 are shown in figEJ This picture clearly 
demonstrates the ultrametric organization of the barriers separating the valleys. 





v(xly=0.008) 


V ..aA.LJ 


L.JLaa.. i 




FIG. 4: Typical shape of the function v{x\y) = — \nf(x) for few values y = 0.04, 0.008, 0.0016. 

As one can see in figlH the distance from the real axis, y, characterizes the number of 
the Cayley tree generation: as closer is the value y to the real axis, as more new generations 
of barriers appear (and as higher are the barriers). The relation of the system of barriers 
shown in figE] with the structure of RSB matrix (|14j) becomes now very straightforward. 
The barriers of smallest scale correspond to the values q\ in the matrix Q(go 5 Qi, Ql), the 
barriers of the next scale correspond to etc... Hence the distance y might be considered 
as the parameter controlling the size of the RSB matrix: as y — > the size of Q tends to 
infinity. Thus, one can conjecture that the number of the replicas, n, is the function of y. 

This relation shall be discussed in more details at length of the next section. The function 
v(x) = — lnf(x), defined on the interval < x < 1, has the properties borrowed from the 
structure of the underlying modular group PSL(2,Z) acting in the half-plane TC(z\y > 0) 
[13, EH]. In particular: 

• The local maxima of the function v(x) are located at the rational points; 
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• The highest barrier on a given interval Ax = [xi,^] is located at a rational point - 
with the lowest denominator q. On a given interval Ax = [xi,X2] there is only one 
such point. The locations of the barriers with the consecutive heights on the interval 
Ax are organized according to the group operation: 

P ±(B P 1 = P 1 ±P 1 (17) 
Qi <?2 qi + <?2 

The figures |U clarifies this statement. The highest barriers on the interval < x < 1 are 
located at the points x = and x\ — 1. Rewriting and 1 correspondingly as ° and \ 
we can find the point Xi of location of the barrier with the next maximal hight. Namely, 
X2 = y © i = i+i = \- Continuing this construction we arrive at the hierarchical structure 
of barriers located at rational points organized in the Farey sequence. According to our 
construction this should replace the standard RSB scheme. The corresponding "modular" 
(MRSB) matrix shall be constructed in the next section. 



C. Explicit organization of barriers on the basis of Dedekind function 

The said above gives the idea of the construction. However to establish the precise 
connection between the standard RSB scheme and the "number-theoretic" MRSB scheme, 
we should construct the MRSB scheme with exponential organization of the Boltzmann 
weights — as in RSB. Recall that the ultrametric distance between two points A and B on 
the Cayley graph is equal to one half of the number of steps of the shortest path connecting 
these points. For example, the ultrametric distances tab (between points A and B) and 
tag (between points A and C) in figOl are as follows: tab = 1 and tac = 3. It is clearly 
seen that the points A and B (as well as the points B and C) are separated by the barrier 
located at the point Xi = |, while the points A and B are separated by the barrier located 
at the point x-i = 1. According to the last of eqs.(|To^. the height U m of the barrier between 
two points separated by the ultrametric distance m in the standard RSB scheme reads 

U* SB = -\nq m = Pm (18) 

The dependence of Boltzmann weights on the ultrametric distance in the MRSB scheme 
should also satisfy (fTHj). This implies that the heights Umrsb of the MRSB-barriers should 
linearly depend on the ultrametric distance between two vertices of the Cayley tree. To 
adjust the heights of the barriers in MRSB and RSB schemes some auxiliary work has to be 
done. The corresponding construction is described in the rest of this section. Begin with the 
standard RSB tree and distribute barriers separating all 2 L boundary points of the Cayley 
tree in the generation L equidistantly in the unit interval [0, 1]. For example, for L = 1 there 
are 2 vertices of the Cayley tree separated by a single barrier located at the point Xi = |. 

The barriers of the second generation L = 2 are placed at the points = \ and x^ = §• 

For L = 3 we have Xg = |, x^ = |, xf^ = |, x^ = | etc. The points x^ = ^ (where j 
is odd and 1 < j < 2 L ) cover uniformly the unit interval [0, 1]. Let us normalize the barrier 
heights as follows. In the finite Cayley tree of L generations there are 2 L " 1 smallest barriers; 
all of them have equal heights. Demand these barriers be of height U(L) = j. Then the 
averaged barriers in some intermediate generation m (1 < m < L) have the height 

L — m + 1 m — 1 , . 

U(m,L) = = 1-— £- (19) 
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Hence, in the first generation m = 1 the height of a single barrier located at the point 
X\ = | is normalized by 1: U(m = 1,L) = 1. This sets our normalization condition. Let 
us note that the choice of the normalization of the barriers might be considered as the 
renormalization of the inverse temperature j3 — > (3L (see (|18|l). The "physical" sense of such 
renormalization consists in the following. When the number of Cayley tree generations L is 
increased, the corresponding landscape acquires more and more small-scale details without 
changing the height of the maximal barrier, which is always normalized by 1. 

Now we proceed in the similar way with the continuous MRSB-structure. Recall that 
the logarithm of the Dedekind function has maxima at the rational points x — |. The set of 
these points {x{\ can be generated recursively as it has been explained above in connection 
with eq. lfTTj) . We start with two "parent" points = % and 1 = ~ of zero's generation. On 
the first step we generate a barrier at the point = on the second step we generate 
two barriers at the points | = 2±1 and | = and so on... The number of recursive steps 
necessary to arrive at a barrier located at a specific point x = | we shall call the generation 
in which this barrier has appeared for the first time. 

Let us fix the maximal generation L. Define {xi}(m) the set of points in the generation 
m < L at which the barriers are located. There are 2 m ~ 1 such points. For example, 
{xi}(m = 1) = {-}, {xi}(m = 2) = {|, §}, {xi}(m = 3) = §, §, §} and so on. Now we 
need to express the heights Umrsb(xi, m, y) of the barriers up to the generation m < L via 
the Dedekind function rj(xi,y). Take the logarithm of the normalized Dedekind function 

h(x,y) = -lnf(x,y) (20) 

as the possible basis for such construction (the function f(x,y) is defined in ©)• Our 
choice is justified by the fact that this function has maxima at the points {xi}(m) and 
demonstrates the ultrametric behavior. Starting from m > 3, the values of the function 
h(xi,m,y) at the points of maxima {xj}(m) become not equal. In fig|S^ we have plotted 
the minimal h min {m,y) = m.m(h(xi,m,y)), the average h avr (m,y) = (h(xi,m,y)) and the 
maximal h max (m,y) = max(h(xi,m,y)) barrier heights in each of three generations m = 
4,6,8. The brackets (...) denote averaging over all barriers in the generation m. The 
function 0.001 y~ l is plotted for comparison in the same figure. 

The figure allows us to conclude that h(m,y) ~ y^ 1 for fixed m. More detailed 
information about the function h avr (m,y) = (—ln(n(xi,y))) can be extracted from the 
figEb, which is similar to the figEt but is extended up to m G [1, 10]. The data shown in 
figlHb is approximated by the natural fit 

havr Q>m,y 

for all points such that h avr > 10. The results of our approximation are shown in the Table 
UK as well as are drawn in figJSjD by lines. 

We have plotted the coefficient a m as a function of m in figEJi. The approximation 
a m ~ 0.144(8) e -°- 847 ( 18 ) m i s shown in the same figure by dashed line. Hence, we expect the 
following relation h(m,y) ~ 0.144(8) e -°- 847 ( 18 ) m y' 1 . 

The analysis of the numerical data for the averaged logarithm of the Dedekind function 
allows us to conjecture the general ansatz 

h avr (m,y) = (-\nf(x l ,y))=d(y)e~ c ^ m y- 1 (21) 



12 




FIG. 5: a) Minimal h m i n (y), average h avr (y) and maximal h T 
function h avr (y) for m = 1, 10. 



• (y) barriers for m = 4,6, 8; b) The 



m 






1 


0.0594(8) 


-1.008(2) 


2 


0.0263(4) 


-1.008(1) 


3 


0.0118(2) 


-1.011(1) 


4 


0.0053(1) 


-1.014(2) 


5 


000230(7) 


-1.021(3) 


6 


0.00084(7) 


-1.04(1) 


7 


0.00031(4) 


-1.06(1) 


8 


0.000118(17) 


-1.08(6) 


9 


0000040(7) 


-1.102(1) 


10 


0000012(3) 


-1.14(2) 



y 


c(y*) 


5.12 x 10~ 5 


0.820(11) 


2.56 x 10~ 5 


0.8057(99) 


1.28 x 10~ 5 


0.7944(86) 


6.4 x 10" 6 


0.7718(29) 


3.2 x 10~ 6 


0.7623(22) 


1.6 x 10" 6 


0.7579(27) 


8 x 10" 7 


0.7556(30) 


4 x 10" 7 


0.7544(32) 


2 x 10" 7 


0.7538(37) 


1 x 10" 7 


0.7535(34) 



TABLE I: a) Numerical results of the fit h avr 
y = y* up to L = 10 generations — see eal231 



a m y bm ; b) Dependence of the parameter c(y) on 



This equation establishes the functional dependence of the average height h avr on the gen- 
eration m (m = 1, 2, ...) for fixed y. Recall that the generation m is defined as the minimal 
number of recursive steps necessary to arrive at a barrier located at a specific point Xi — |. 
Let us fix now the maximal generation L. Our desire is to obtain the normalized barrier for 
MRSB scheme like it has been done for RSB (see eq. (|19j) ). Taking into account that ()21j) is 
valid for any m > 1, we may extract m from (}2l"| and write it as follows 



m 



1 + 



h avr (m,y) 
c{y) "~ h avr (m = l,y) 



In 



(22) 



Substituting (|2*2*|) into (|T9*|) we get the following formal expression for the height of the 
barriers in the MRSB scheme 



U M RSB(x m ,y,L) = 1 



m 



1 



1 + 



h aV r(m,y) 
c(y)L"^ h avr (l,y) 



In 



(23) 
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FIG. 6: a) Approximation of the coefficient a m ; b) Ultrametric barriers UMRSBi^iV) = 
i !„ K m ,y)) 

Now we should explain how it is possible to use the formal equation ()23|) in practical com- 
putations. Let us fix the maximal generation L and define the value y* from the equation 

h avr (y*,L) = l (24) 

The value y* has the following sense. For all y < y* the value h avr (m,y) does not depend 
on m (m < L) because all barriers in the generation m are already presented and hence the 
average hight h avr can not be changed for any y < y*. Therefore we can postulate, that for 
all y < y* the height of the barrier is described by the function 

UmbsbM = 1 + In ^ f }^ V 2 (25) 

c(y*)L ln/(|,y*) 

where Xi is the rational point of the barrier location and we have used (j2*Uj) . Special attention 
deserves the discussion of the dependence c(y*). As it is shown in the Tabled, for sufficiently 
small y* the value of c tends to the constant value. Thus the function UMRSB^i) defines 
the heights of ultrametric barriers in MRSB scheme fully consistent with the heights of the 
RSB-barriers for the Cayley tree. It should be noted that for intermediate values of y some 
barriers of large generations are negative due to the cut (J24|) . however having m, we can 
always find y* such that all barriers are positive. For example, for y < y* = 10~ 6 all barriers 
up to m = 10 generations are positive. 

Our tedious but simple construction is clearly explained in the figure [7[ where we have 
plotted the heights of barriers UMRSB^Xi) at the rational points x« up to the generation 
L = 10 for y* = 10~ 6 . Once the maximal generation L is fixed and all barriers for all 
generations (up to the maximal one) are shown in the picture, then this picture is unchanged 
for all y such that y < y*. According to our definition, these barriers are normalized in such 
a way that the height of the largest barrier is equal to 1 . The averaged height of the barriers 
of second (third, forth, etc..) generation is equal to 0.9 (0.8, 0.7, etc..) - compare to ()19j). 
Let us repeat that in fig|7|we have displayed the barriers corresponding to generations m up 
to m — L. We can easily relax this condition and work with the set of all possible barriers 
up to some fixed minimal value y m in without paying attention how many generations are 
counted. 
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FIG. 7: Ultrametric barriers Umrsb( x i) at the rational points X{ for y = 10~ 6 

D. Description of dynamical models 

Our aim is to investigate numerically the probability distribution of a random walk on 
the boundary of an ultrametric space. This problem has been a subject of the original work 
[l5| where the dynamics in the ultrametric space was considered in the frameworks of the 
discrete model of jumps on the boundary of 3-branching Cayley tree. The first application 
of p-adic analysis for the computation of the probability distribution of the diffusion at the 
boundary of the ultrametric tree has been successfully realized in |6(. 

Below we reconsider the diffusion problem in the continuous metric space, where the 
ultrametric organization of the barriers is due to the modular structure of the Dedekind 
function. The advantages of such consideration we shall discuss in the Conclusion. Let us 
specify the models under consideration. 



1. I. The standard RSB model 

In this model we consider the random jumps at the boundary of 3-branching Cayley 
tree — see figOl The transition probabilities between the Cayley tree vertices are encoded in 
the matrix ()14|) . As explained above, we fix the total number of generations L and set the 
highest barrier (in the first generation) always equal to 1. Then the heights of the barriers 
of generation m are given by (fT^j) . while the height of the smallest barriers in the generation 
L is equal to Increasing the number of generations L, we increase the resolution of our 
model without changing the value of the maximal barrier. We use this problem of diffusion 



15 



on the boundary of the standard 3-branching Cayley tree as a testing area for diffusion in 
the continuous tree-like ultrametric space. 

Let Pisf(k, ko) be the partition function of the random walk starting at the point ho and 
ending after N steps at the point k, where k labels the vertices of the fc's generation of 
the Cayley tree and both points k^ and k belong to the m's generation of the tree. The 
corresponding diffusion is governed by the recursion relation 



P^(k)=J2Q RSB (k,k')P^(k') 

k'=0 

where Q RSB (k,k') is the 2 L x 2 L matrix whose elements are the Boltzmann weights asso- 
ciated with jumps from some point k\ to the point k 2 ({£4, k 2 } G [1,2 ]) at the boundary 
of 3-branching Cayley tree cut at the generation L. In order to take into account the ul- 
trametricity of the Cayley tree, the matrix Q RSB (k, k') should have the RSB structure and 
hence should coincide with the matrix Q in (|14j) . 



2. II. The quasi- continuous MRSB model 



This model is defined as follows. We fix the maximal generation L and consider the heights 
of the potential JJ M (see (I25|0 corresponding to generations less or equal L. For example, 
in figlHt we have plotted 15 = 2 L — 1 = n barriers up to L = 4 generation. The rational 
points Xi split the interval [0,1] into 16 subintervals [xi_i,Xi], i = 1,2, ...,15 (including 
the subintervals [0,a;i] and [#15, 1]). The random walker can jump form one interval to any 
other. Let us enumerate the intervals by the variable A; = 1,2, 3, ...,16. So, k — 1 denotes the 
interval [0, x\], k = 2 is the interval [xi,^], etc. Let UMRSB^k, k') be the maximal barrier 
between the intervals k and k' defined by (J2HJ). By construction U MRSB (k, k) = 0. According 
to our rules if x G [xk-i,Xk\ and x' G [x^'-i, x^] then U MRSB (x,x') = U MRSB (k,k'). Thus 
we can write the recursive relation for the continuous distribution function Wn{x). This 
function is defined for x G R, (0 < x < 1), where iV enumerates the time moments. We 
claim for this quasi-continuous MRSB case the following form of the master equation which 
should replace eq. (f2H . 

pi 







W N=0 (x) = S : 



We can transform (J27j) making it more similar to (JUJ). Namely, define P^ RSB (k) - the 
probability that the walker belongs to the A;-th interval on the iV-th jump (time step) 

Px IRSB (k) = W N (x)dx (28) 

The normalization of the probability defines the sum over all intervals Xwb=i PN RSB {k) — 
Jo Wn(x)cIx = 1. Integrating (j2T|) over the interval [0, 1] and replacing the integrals over dx 
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FIG. 8: a) Location of barriers up to 4 generations: a) Model I: the values v i correspond to the 
intervals between maxima; b) Model II: the values Vi correspond to the discrete points. 



and dx' by sums over k and k', we get 

- N+ i(x)dx 



k=i Jo 



1 r l 



JO 



x k fX 



e -PU M RSB(x,x') WN ( x ') dx 'dx 
= J2J2I l K e-^ U MBSB{x,x') w MRSB^ )dxldx ^ 

k=l k'=l J x k-1 J x ' k '-1 
n n 

= rr i ke -f3U M RSB(k,k') P MRSB 
fe=i k'=i 

where 4 is the length of the k-th interval 4 = f Xk dx = Xk — %k-i- Here we have used (j2lfl) 

to express P^f B (k) and P^ RSB {k') via W^f B (x) and W$ RSB (x'). Now we can separate 
equations for different k and write 

n 

P^ + T(k) = W. k>)P*? RSB (k') (30) 



k'=l 



where 



Q(k, k') = l k e 



-PU MRSB (k,k') 



(31) 



is the probability to jump from the interval k' to the interval k. This probability is 
the product of the measure of the appropriate interval Ik and the Boltzmann weight 

e~@U (k,k) Define a vector P^ RSB whose elements P^ IRSB (k) are the probabili- 
ties for a walker to reach the interval k on the iV's jump. The vector p^ IRSB is the initial 
probability distribution on the intervals. Introduce a matrix Q with the elements Q(k, k'), 
and a matrix \J MRSB with the elements U M R B (k,k'). The elements of these two matrices 
are connected by the relation ()31j) . The probability for the walker to reach the interval k at 
the time moment N is the fc's element of the vector PfJ RSB (k) 



,MRSB 

N 



(k) = Q J X 



N-cMRSB 



(32) 
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Let us remind, that the barrier U MRSB (k,k') between the points k and k! is the highest 
barrier on the interval [k, k'\. Such barriers are located at the rational points {k m }. This 
defines the ultrametric structure of the matrix \J MRSB and, therefore, of Q. By construction 
the matrix \J MRSB reflects the ultrametric hierarchy of barriers appearing in the standard 
RSB-matrix (see eq. ()14j) ). For example, the structure of the "modular" MRSB matrix 
jjmrsb defying the heights of the barriers between all the intervals along the x-axis for 
L = 8 is displayed below: 



U 



MRSB 



Uo U(i) 
U o 


u(i) u(l) 
u(l) u(l) 


u(l) u(l) u(i) u% 

U(l) U® 17(1) 17(1) 

u % u % 

17(1) ^(1) i7(|) 1/(1) 


u(l) u(l) 


Uo U(l) 

u(l) u 


U® U% 17(1) i7(§) 

u (l) u% u (|) [/(I) 

^(1) u% 
u% £/(§) 


U(l) U 


u(l) u(l) 
u(l) u(l) 


u(l) u(D 
u(l) u(l) 


Uo U(i) 



3. III. The discretized MRSB model 



In this model the location of the barriers is the same as in the model II, but now we split 
the total interval [0, 1] into n = 2 L sub intervals of the equal length - each. The midpoint 
of the k-th interval is the point Xk = — 1- + -. We consider the distance between intervals 
k and k' as the distance between the midpoints: Uin(k,k') = U M (xk,Xk>). Namely, we 
consider jumps between n midpoints k, (k = 1, 2, n) with coordinates Xk = — ^ + -■ The 
probability to jump from a point k' to a point k is defined by the relation 

Qiii^ k ' ) = l e -PU MRSB i^y) (34) 
L 

Each element P^ n (k) of the of the probability vector P corresponds to a point. In this 
model all barriers Xi = — lie between the midpoints — see figJHjD. 

Comparing the models II and III it is easy to understand, that the number n of points 
Xk in the interval between barriers [xk, Xk+i] plays role of measure of that interval. 
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Below we compare the numerical solutions of the models I, II and III. To be more specific, 
we consider simultaneously the iV-step random walk on the Cayley tree (the model I) and 
on its continuous analog (the models II and III). First of all, we calculate the elements 
[Q N ](i,j) of the matrix for all three models. Then we compute: 

- The probability Pq to stay at the initial point after N random jumps 

Po^E^KM) (35) 

i=l 

- The average ultrametric distance U between the initial and the final points of the 
iV-step random walk: 

n n 

^EEi u ]^')[Q w ]^') ( 36 ) 

i=l j=l 

where [U](i,j) is the element of the matrix U defining the ultrametric distances be- 
tween all available points. The size of the matrices U and Q are n = 2 L , where L is 
the total number of generations. 

E. Numerical results for the models I, II, III 

The following three cases are studied numerically: 

1. We fix the maximal generation, L, and the number of jumps, N, and plot P an d U as 
the functions of the inverse temperature (3 - see figsJHKjb for the following numerical values 
of L and N: L = 10, N = 400. 




P P 
FIG. 9: a) Pq(0); b) U(Jf). Numerical values of the parameters are: L = 10, N = 400. 

2. We fix the maximal generation, L, and the inverse temperature, (3, and plot Pq and U as 
the functions of the number of steps iV - see figsllOb.b for the numerical values of L and (3: 
L = 10, P = 4. 

3. We fix number of jumps, N, and the inverse temperature, (3, and plot Pq and U as the 
functions of the maximal generation L - see figsJTTk.b for the following numerical values of 
N and (3: N = 400, [3 = 1. 
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FIG. 10: a) Pq(N); b) U(N). Numerical values of the parameters are L = 10, /3 = 4. 
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FIG. 11: a) Pq(L); b) U(L). Numerical values of the parameters are N = 400, = 1. 



III. DISCUSSION 



The numerical solution of the master equations ()26|) and ()27|) shown in figsl9TITTl allow 
us to conclude that we arrive to the same results for the return probability Po an d the 
mean height of the barriers U for random walks in the ultrametric spaces independent on 
which model we consider: either the ultrametric structure of the barriers is defined on the 
boundary of the Cayley tree (case I), or the heights of these barriers are given by the properly 
normalized Dedekind function (cases II and III). 

Our results are obtained under the condition that the heights of the ultrametric barriers 
in all models (I, II and III) have linear dependence on the number of Cayley tree generation. 
Recall that the construction of our ultrametric potential (|25|) on the basis of Dedekind 
modular function is conditioned only by this demand. The current choice of the potential 
is stipulated basically for demonstrational aims. Namely, we have shown that it is possible 
to adjust the ultrametric structure of the continuous metric space 7i to the structure of the 
standard Cayley tree in such a way that we can mimic the main statistical properties of the 
random walk at the boundary of the Cayley tree by the diffusion in Ti. 
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The key ingredient of our construction is contained in the replacement of the RSB scheme 
(JHJ) by the new MRSB scheme Despite yet the corresponding master equation (|27jl 

allows only numerical treatment, the proposed construction, to our opinion, has some ad- 
vantages with respect to the RSB scheme resulting from: i) the self-duality (see eq. ()13p ). 
and ii) the continuity of MRSB. 

We permit ourselves to touch below both these properties (i) and (ii) in connection with 
the possible simplification of the solution of the continuous analog of the master equation (J29J) 
(see part IIIlAI below), and with the speculation about the geometry of n — > limit (see part 
IIII Bl below). We clearly realize that the last subject is one of the mostly intriguing question 
is the statistical physics of disordered systems and our conjecture might be criticized from 
different point of view, however we believe that our geometric construction might stimulate 
the readers to fill the standard mode of thinking about the replicas by some fresh geometric 
sense. 



A. The continuous MRSB model 



Taking advantage of our construction of ultrametric system of barriers described by the 
potential UMRSB{.x m ) (see ()25|)). we can replace the master equation (pjj) by its analog based 
on the properties of the Dedekind function. The distribution function W^{x) is described 
by ()27|) where x m is the rational point on the interval [x', x] with the smallest denominator 
(there is only one such point). As we have already seen, the highest barrier UMRSB^m) 
on the interval [x',x] is located at the point x m . The number of barriers and their heights 
are controlled simultaneously by the parameter y. As y — > 0, more and more new barriers 
appear. Gradually they become higher and narrower preserving the ultrametric structure. 
In this section we are interested in the limit y — > 0, where the hierarchical structure of 
barriers is highly developed — see, for example figd Simplification of eq.([27|) is allowed by 
using eq.(JT3J). Making in eq. (}2Tj) the substitution which maps y to ^ we get: 



tt i t*\ - 1 j 1 i ln f( x m,y) 

U M RSB{Xm, V, L ) = 1 + In 



cL*(y) Infdy) 



cL*(y) q 



In - (37) 



where L*(y) is determined again using the relation (|23J)- The equation (pTTj) is valid only for 
{p, q,s,r} G Z such that ps — qr = 1 (where r might be an arbitrary integer). The point 

| j and its "dual image" x' m = ||| lie in the unit interval [0, 1]. Substituting (jHZj) 



■Em 

into P7jl . we have: 



rl /n\W(^*(!/)) 

W N+1 (x) = J dx'e-Pl-J W N (x') 



(38) 



W N=0 (x) = 8, 



X,Xg 



Eq.()38J) has rich number theoretic structure. In fact, the information of the position of the 
highest barrier on the interval [x', x] is hidden in the integer q. Let us recall that the value 
of q should satisfy the following condition: q is the lowest denominator of the rational point 
x m = 2 x m G As one sees, the value L*(y) governs the amplitude of the potential, 

while q controls the maximal "resolution" and, hence, the total number of barriers. 
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B. Speculations about continuous number of Cayley tree generations L and replica 

n — > limit 

Despite of many advantages of the RSB scheme, it contains the unavoidable in the replica 
formalism mystery of the n — > limit, where n = 2 L is the number of replicas — see fig ED 




54 7 38572 7 5 8 3 7 4 5 



FIG. 12: The 3-branching Cayley tree. The barriers are located at the rational points ordered in 
the Farey sequence. The the replica limit n — * is interpreted as the absence of states (vertices of 
the tree) when y — ► oo. 

One of the merits of our construction based on the Dedekind ^-function consists in the 
possibility to change continuously the number of ultrametric barriers and their heights by 
varying the parameter y, and hence to change continuously the number of states n of the 
RSB matrix. Let us remind that the standard 3-branching Cayley tree being isometrically 
embedded in the space 7i = {z\ Imz > 0} becomes a "continuous" structure: one can define 
the smoothed neighborhood of maxima of the function f(z) in the half-plane y = Imz > 0— 
see fig[D For instance, our construction allows us to increase continuously the number of 
states n from, say, n = 1 to n = 2 changing correspondingly the parameter y in the interval 
[yn=i,Vn=2\, where y n=1 = and y n=2 = ^f. 

Our geometric construction based on the isometric embedding of the Cayley tree into the 
complex plane permits us to consider also the opposite limit — the case when n is formally 
less than 1. This case means the absence of any states in the RSB matrix (n = 0). As 
one sees from figsfflandEl the function f(z) has no longer local maxima above the value 
y = Im z > As soon as we have identified the local maxima of the function f(z) with the 
vertices (states) of the Cayley tree, we conclude that the absence of any states of RSB matrix 
means the absence of local maxima of the function f(z). Thus, the limit n^Owe interpret 
as y — > oo. This claim is based on the behavior of the function f(z) in the upper half-plane. 
The correspondence of the limits n — ► and y — > oo is schematically illustrated in the figfT^l 
The geometrical interpretation of number of replica states in the region < n < 1 makes 
the physical statement about the "analytic continuation n 0" more formal. Actually, let 
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us recall the definition of the hyperbolic distance £ in Ti(z\Im.z > 0) between two points 
z = (x ,y ) and z = {x,y): 

coshC=l+ iX - Xo)2 + iy - yo)2 (39) 

yoy 

As we have seen, the hyperbolic distance £ is the continuous analog of the number of Cayley 
tree generations, L, what permits us to define the number of replica states n as 

n = e constC (40) 

As y — > 0, we can replace cosh£ with the leading exponential term. This way we arrive at 
the standard relation e c = — . Hence, 



yo 



e , \ = ( x - x o) + Vl (ai) 

2 T voy 



neglecting the second term in the left-hand side of ([41)1. we get 

£ = -\ny + d (y -> 0) 

where Cl = In 2((a - 3!o)2+yg) . 
1 yo 



However if y — > oo, eq.()39p gives 



— = ]L (42) 
2 y 1 ; 



where according to the physical condition n —>■ and eq. ()40j) . we have to take another branch 
of the function arccosh (...), corresponding to the negative values of £. This leads to the 
following definition of £: 

£ = -\ny + c 2 (y^oo) (43) 

where c 2 — In y. As y — > oo, the hyperbolic distance £, defined according to (jjHj) tends 
to — oo. Equation ()43j) is consistent with the definition (}4T)j) of number of replicas n in the 
region < n < 1. 

It would be very desirable to check how our conjecture works in the models possessing 
the RSB symmetry of the order parameter. We expect that our construction might be useful 
for the investigation of aging phenomena considered from the point of view of the diffusion 
in the whole ultrametric space Ti.(z\lmz > 0). 
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